Search for annual modulation in low-energy CDMS-II data 
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We report limits on annual modulation of the low-energy event rate from the Cryogenic Dark 
Matter Search (CDMS II) experiment at the Soudan Underground Laboratory. Such a modulation 
could be produced by interactions from Weakly Interacting Massive Particles (WIMPs) with masses 
~10 GeV/c^. We find no evidence for annual modulation in the event rate of veto-anticoincident 
single-detector interactions consistent with nuclear recoils, and constrain the magnitude of any 
modulation to <0. 06 event [keVnr kgday]"^ in the 5-11.9keVnr energy range at the 99% confidence 
level. These results disfavor an explanation for the reported modulation in the 1.2-3.2 keVco energy 
range in CoGeNT in terms of nuclear recoils resulting from elastic scattering of WIMPs at >98% 
confidence. For events consistent with electron recoils, no significant modulation is observed for 
either single- or multiple-detector interactions in the 3.0-7.4 keVco range. 

PACS numbers: 14.80.Ly, 95.35. +d, 95.30.Cq, 95.30.-k, 85.25.Oj, 29.40.Wk 
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Astrophysical observations indicate that the vast ma- 
jority of matter in the universe consists of non-baryonic, 
non- luminous dark matter [T]. Weakly Interacting Mas- 
sive Particles (WIMPs) are leading candidates for this 
dark matter since they would be thermally produced in 
the early universe in the correct abundance to account 
for the observed relic density [H |2]. Possible experi- 
mental signals from DAMA/LIBRA [3], CoGeNT gj and 
CRESST-II [S] have led to significant recent interest in 
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The WIMP scattering rate is expected to annually 
modulate due to the relative motion of the earth through 
the local dark- matter halo [9]. The presence of an an- 
nually modulating component in the observed interac- 
tion rate can identify a WIMP signal in the presence 
of significant unmodulated backgrounds. This modula- 
tion signature is especially useful for WIMPs with masses 
~10GeV/c^ which would primarily produce recoils with 
energies just above the detection threshold, where the 
rejection of backgrounds that can mimic a WIMP sig- 
nal is less powerful. Both the DAMA/LIBRA [3] and 



CoGeNT [1] experiments claim evidence for such a mod- 
ulating signal in their data. 

The CDMS II experiment consists of an array of cryo- 
genic germanium and silicon detectors which measure 
both the ionization and athermal phonon energy de- 
posited by each particle interaction [101 HI]- The ra- 
tio of the ionization to phonon energy allows discrimi- 
nation of the expected WIMP nuclear-recoil signal from 
more prevalent electron-recoil backgrounds on an event- 
by-event basis. At the <10keV recoil energies expected 
from WIMPs with masses <10GeV/c^, electron-recoil 
rejection is less effective, because the ionization signal 
is comparable to readout noise. Even with this limited 
rejection, the observed interaction rate in the CDMS II 
germanium detectors has been shown in a previous analy- 
sis [12 to disfavor an explanation for the DAMA/LIBRA 
and CoGeNT signals in terms of spin-independent elastic 
scattering of low-mass WIMPs. 

If, as recently suggested, only a small fraction of the 
low-energy excess events in CoGeNT are due to WIMPs, 
then constraints from CDMS II may be avoided [51 [13] . 
In this case, if CoGeNT's annual modulation is due to 
WIMPs, the fractional variation is several times larger 
than expected for a "standard" halo with a Maxwellian 
velocity distribution [6l [14] . In addition, the energy spec- 
trum of the modulation extends to higher energies than 
expected for a standard halo [Tl lTiHTB] . Such large modu- 
lation fractions and hard spectra might be possible if the 
halo exhibits (non-standard) local substructure [5] 13 [H] . 
To test such a scenario, this Letter searches for a cor- 
responding annual modulation in the CDMS II germa- 
nium data. This analysis does not cover the full energy 
range of the CoGeNT modulation, restricting itself to en- 
ergies above 5 keVnr (which, due to quenching [TT] [TH] , 
corresponds to 1.2keVee in the standard CoGeNT en- 
ergy scale [19]). Because germanium serves as the target 
material for both CDMS II and CoGeNT, these results 
provide a check of whether the reported modulating sig- 
nal is due to WIMPs that is less model-dependent than 
recent results from XENONIO ^ and XENONIOO [IJ. 

The data analyzed here were collected over nearly two 
annual cycles, from October 2006 to September 2008, 
using all 30 Z-sensitive Ionization and Phonon (ZIP) 
detectors installed at the Soudan Underground Labora- 
tory [ini m] . Data-quality and detector selection criteria 
are identical to the previous analysis of the low-energy 
CDMS II nuclear-recoil spectrum in [12]. Only the 8 
germanium detectors with the lowest trigger thresholds 
were used to search for WIMP interactions, while all 30 
detectors were used to veto events with interactions in 
multiple detectors. 

Following [12], the nuclear-recoil energy scale is based 
on the phonon measurement, which is corrected by ~20% 
to take into account the fraction of the total phonon 
signal arising from the Neganov-Luke phonons [22] gen- 
erated by the charge-carrier drift across the detectors. 



The Neganov-Luke phonon contribution for nuclear re- 
coils in ^^^Cf calibration data was directly measured for 
the recoil-energy range of this analysis. As in [12], the 
phonon energy scale for electron recoils was conserva- 
tively calibrated, ensuring to the 90% C.L. that the 1.3 
and 10.4 keV activation-line energies were not underesti- 
mated. 

The maximum energy considered in this analysis was 
11.9keVnr, matching the highest energy observable by 
the CoGeNT "LG" (low gain) channel [H]. Because 
time-dependent variations in the CDMS trigger thresh- 
olds could mimic or hide a modulation in the event rate, 
the energy threshold for this analysis was conservatively 
chosen to be SkeVnr, high enough that events are trig- 
gered with essentially perfect efficiency. To avoid bias, 
trigger efficiency was measured throughout the exposure, 
using events for which at least one other detector trig- 
gered. For the 5-11.9 keVnr energy range, combining 
all detectors and all time bins yielded 4350 events in 
this unbiased sample, only 3 of which failed to trigger. 
These missed triggers were each in a different detector 
and were uniformly spaced throughout the considered en- 
ergy range. 

Because CDMS II uses a phonon-based energy scale 
(at these low energies), and CoGeNT uses an ionization- 
based energy scale, quenching causes the two experi- 
ments to exhibit different mappings between energies 
assuming nuclear recoils (such as the energy range of 
this analysis, 5.0-11.9 keVnr) and energies assuming elec- 
tron recoils. For electron recoils with the same total 
phonon signal in the CDMS II experiment, the equiv- 
alent recoil-energy interval is 3.0-7.4 keVoo, due to the 
larger Neganov-Luke phonon contribution. Analogously, 
for electron recoils with the same total ionization signal 
in the CoGeNT experiment, the equivalent recoil-energy 
interval is 1.2-3.2 keVoc, where we apply CoGeNT's mea- 
sured ionization yield for nuclear recoils [121 US] • 

Detector stability was monitored throughout data tak- 
ing with quality cuts, removing periods of abnormal de- 
tector performance jllj . For consistency with previous 
work, we followed [12] and removed data taken during 
the 20 days following exposure of the detectors to a neu- 
tron calibration source. After removing these time peri- 
ods, a total of 241 kg days raw exposure were considered, 
as in [12]. To allow checks for stability to be applied to 
multiple-scatter events, an additional cut was introduced 
to eliminate electronics "glitch" events, for which phonon 
pulses were detected above threshold in more than 15 de- 
tectors simultaneously. 

Events inconsistent with WIMP interactions were re- 
jected. Since modulation of data-selection cut efficiencies 
could mimic or hide a modulation in the event rate, selec- 
tion criteria were designed to have constant acceptance 
with time, and any residual modulation in the cut effi- 
ciencies was constrained using events sampled through- 
out the data taking period. Since WIMPs have a negligi- 



bly small probability of interacting more than once in the 
apparatus, events with energy deposited in more than a 
single detector ( "singles cut" ) or in the active scintillator 
veto ("veto cut") were removed. The glitch cut, veto cut 
and singles cut have a combined efficiency >97%, with 
negligible time-dependent variation. 

Events were further required to have ionization sig- 
nals consistent with noise in the outer charge electrode 
of the detector ( "Q- inner cut" ) . To search for a nuclear- 
recoil signal of WIMP origin, the ionization energy was 
required to be within ±2cr of the mean of the energy- 
dependent nuclear-recoil distribution from calibration 
data ( "nuclear-recoil cut" ) in the main analysis described 
below. This ionization-based selection increases the sen- 
sitivity of this analysis to a modulating signal relative 
to the more restrictive ionization-based selection used 
in [12], provided that backgrounds do not modulate. To 
explore different physics or instrumental origins of a po- 
tential signal, we also applied our modulation analysis to 
two additional event samples consisting of either single- 
scatter or multiple-scatter events with no ionization- 
based nuclear- recoil cut. The quality, glitch, veto and 
Q-inner cuts were always maintained. 

These various cuts result in an efficiency e (i, E, d) 
that depends on the time t, the deposited energy E and 
the detector d. With notations following the descrip- 
tions above, the total efficiency for our primary "WIMP- 
candidate" sample can be written as 

£ [Z^ rj ^ Qjj Cgiitch^triggcr^singlcs v^/ ^vcto ^ V^y 

£Qinnor (^i E, d) ENuclcarRccoil {t, E, d) , 

where we have explicitly identified the dependence on 
time t, energy E and detector d for each of the cuts. 
For the event samples that remove the nuclear-recoil and 
singles cuts, the corresponding efficiencies, Esingics and 
ENucioar Recoil, are ignored. 

In the lower portion (5-7.3 keVnr) of the energy range 
considered, the rate of nuclear-recoil candidate events 
measured in this analysis is 0.28±0.03 [keVm- kg day]~^, 
while the maximum-likelihood estimate for the CoGeNT 
modulation amplitude is 0.35 [keVnr kg day]~^. The cor- 
responding numbers for the entire energy interval con- 
sidered are 0.15±0.01 [keV„r kg day]"! for CDMS and 
0.16 [keVnr kg day]^^ for CoGeNT. In both cases, a mod- 
ulation of the magnitude observed by CoGeNT would 
require a modulation fraction in CDMS of ~100%. 

We test whether the Q-inner and nuclear-recoil cut ef- 
ficiencies are sufficiently constant using calibration data 
collected throughout the time period used in the analysis. 
For each cut, for a given time interval 7, and detector d, 
we measure P-yd events passing the cut and F-yd failing. 
Note that the time intervals of the efficiency data are not 
coincident with the low-background time intervals, but 
are suitably distributed over the whole data-taking pe- 
riod. We then maximize the likelihood appropriate for a 
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The efficiencies e^d are written as 
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where w = 27r/365.24day~^. For the chosen cut and 
energy interval we fit for the detector-dependent un- 
modulated efficiency e^, the detector- independent rela- 
tive modulation amplitude A, and the phase (p (measured 
from Jan. 1st), while requiring the efficiency to be <1. 
We generate 10^ artificial realizations of the model un- 
der consideration, and determine the confidence regions 
for the modulation and phase using the Feldman-Cousins 
method |24j . This analysis indicates that the maximum 
efficiency modulation allowed by our experimental mea- 
surement of the nuclear-recoil cut efficiency is 1.2% at 
the 90% C.L. In the case of the Q-inner cut, this upper 
limit is 2.3%. 

In order to estimate the modulation in the observed 
event rate, we bin events into 16 time intervals, labeled 
by /?, of ~25 days each. We denote the center of each 
time bin as tp and its width as Ai^. The number of 
events observed in the time interval /3 in detector d is 
npd- We construct a likelihood using the expected Pois- 
son distribution for the n^d 
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where factorial terms have been omitted for convenience. 
In this equation, /i^^ is the expected number of events 

M;3d = {r^ + M cos [w {tfi - 0)]} mdefjdffSd^tfjAE, (5) 

where F^ is the unmodulated rate in detector d, M is 
the modulation amplitude, nid is the mass of detector d, 
epd is the appropriate efficiency using Eq. [ll fjjd is the 
live-time fraction appropriate for the detector and time 
interval, Ai^ is the time interval width, and /S.E is the 
energy interval width. 

Since the efficiency modulation allowed by the fits to 
Eq. [3]is much smaller than the physics effect we are test- 
ing, we need not add an additional term ^gff (^/3d) in the 
likelihood, which would take into account such uncertain- 
ties in Ec^. [31 

Figure [T] shows residual rates for WIMP candidate 
events, after subtracting the best-fit constant rates F^ 
(found with modulated rate M fixed at 0). Using a 
Feldman-Cousins approach, we test modulation models 
[M , 0] on the WIMP candidates, which consist of all 
events satisfying the data-selection cuts described above. 
Figure [2] shows that our observed WIMP-candidate event 
rate is consistent with a constant value. All modulated 
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FIG. 1. (color online) The rate of CDMS II nuclear-recoil 
band events is shown for the 5.0-11.9 keVnr interval (dark 
blue), after subtracting the best-fit unmodulated rate, Fd, 
for each detector. The horizontal bars represent the time 
bin extents, the vertical bars show ±1(t statistical uncertain- 
ties (note that one CDMS II time bin is of extremely short 
duration). The CoGeNT rates (assuming a nuclear-recoil en- 
ergy scale) and maximum-likelihood modulation model in this 
energy range (light orange) are shown for comparison. The 
GDMS exposure starts in late 2007, while the CoGeNT expo- 
sure starts in late 2009. 



rates in this energy range with amplitudes greater than 
0.06 [keVnr kg day]-i are excluded at the 99% C.L. 

For comparison, a similar analysis was carried out us- 
ing the publicly available CoGeNT data (TH] . Our analy- 
sis of CoGeNT data is consistent with previously pub- 
lished analyses O [7j [14]. Figure [s] shows the modu- 
lated spectrum of both CDMS II and CoGeNT, assum- 
ing the phase (106 days) which best fits the CoGeNT 
data over the full CoGeNT energy range. Compatibil- 
ity between the annual modulation signal of CoGeNT 
and the absence of a significant signal in CDMS is de- 
termined by a likelihood-ratio test, which involves cal- 
culating A = Cq/Ci, where Co is the combined max- 
imum hkelihood of the CoGeNT and CDMS data as- 
suming both arise from the same simultaneous best-fit 
values of M and (/>, while £i is the product of the maxi- 
mum likelihoods when the best-fit values are determined 
for each dataset individually. The probability distribu- 
tion function of —2 In A was mapped using simulation, 
and agreed with the x^ distribution with two degrees 
of freedom, as expected in the asymptotic limit of large 
statistics and away from physical boundaries. The simu- 
lation found only 82 of the 5x 10'^ trials had a likelihood 
ratio more extreme than was observed for the two ex- 
periments, confirming the asymptotic limit computation 
which indicated 98.3% C.L. incompatibility between the 
annual-modulation signals of CoGeNT and CDMS for the 
5.0-11.9 keVnr interval. 

We extend this analysis by applying the same method 
to CDMS II single-scatter and multiple-scatter events 
without applying the ionization-based nuclear-recoil cut. 
These samples are both dominated by electron recoils. 
Figure |4] shows the confidence intervals for the allowed 
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FIG. 2. (color online) Allowed regions for annual modulation 
of CoGeNT (light orange) and the CDMS II nuclear-recoil 
sample (dark blue), for the 5.0-11.9 keVnr interval. In this 
and the following polar plot, a phase of corresponds to Jan- 
uary 1st, the phase of a modulation signal predicted by generic 
halo models (152.5 days) is highlighted by a dashed line, and 
68% (thickest), 95%, and 99% (thinnest) C.L. contours are 
shown. 
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FIG. 3. (color online) Amplitude of modulation vs. energy, 
showing maximum-likelihood fits for both CoGeNT (light or- 
ange circles, 68% confidence interval shown with vertical line) 
and CDMS nuclear-recoil singles (dark blue rectangles, 68% 
confidence interval given by rectangle height). The phase that 
best fits CoGeNT over all energies (106 days) was chosen for 
this representation. The upper horizontal scale shows the 
electron-recoil-equivalent energy scale for CoGeNT events. 
The 5-11.9keVnr energy range over which this analysis over- 
laps with the low-energy channel of CoGeNT has been divided 
into 3 (CDMS) and 6 (CoGeNT) equaLsized bins. 
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FIG. 4. (color online) Confidence limits on the amplitude and 
phase of annual modulation for two electron-recoil-dominated 
data samples: multiple scatters (light blue) and single scatters 
(dark red), as defined in the text for the interval 3.0-7.4 keVcc- 
These events are of the same total phonon energy (recoil -f 
Neganov-Luke) as the nuclear-recoil band events of the main 
modulation analysis shown in Fig. [2l of 5.0-11.9 keVnr. 



modulation amplitudes and phases for these two samples, 
both of which are consistent with no modulation. For the 
energy range chosen for this analysis, there is not sig- 
nificant overlap with the corresponding CoGeNT energy 
range under the hypothesis of an electron-recoil modula- 
tion. Our minimum electron-equivalent energy is 3 keVoo 
compared to a 3.2keVcc maximum energy for the Co- 
GeNT low-energy channel. Consequently, this analysis 
cannot exclude the possibility of the modulation observed 
by CoGeNT being the result of electron recoils. The ab- 
sence of modulation in the single-scatter and multiple- 
scatter events indicates the absence of strong systematic 
effects in our data. 
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Supplementary Material 

The following pages provide relevant additional information. 



Energy Scales 

The total energy deposited by a particle interaction is the "recoil energy," En. The majority of the recoil energy is 
deposited directly as athermal phonons, while a fraction of the energy goes into producing electron-hole pairs above 
the band gap. After the charge carriers are drifted across the detector and relax to the Fermi level, the remaining 
recoil energy is deposited in the phonon system. 

On average one electron-hole pair is produced for every e = 3.0 eV of recoil energy for an electron recoil in Ge P51l26| . 
The "ionization energy" Eq is defined for convenience as the recoil energy inferred from the detected number of charge 
pairs, Nq, assuming that the event is an electron recoil: 

Eq = NqX e. (6) 

The ionization energy scale is calibrated in situ using electron-recoil lines of known energy, and is reported in units 
of "keVoo," which gives the recoil energy in keV for an electron recoil producing the same ionization signal. The 
ionization yield Y = Eg/Efi, so on average F = 1 for electron recoils. Nuclear recoils produce fewer charge pairs, 
and hence less ionization energy, Eq, than electron recoils of the same recoil energy. The ionization yield Y for 
nuclear-recoil events depends on recoil energy, with Y « 0.2-0.25 in Ge for Efi « 5-10 keV [T71ITS] . 

An additional phonon population is produced by drifting the charge carriers across the detector, with the total 
energy deposited equal to the work done by the electric field. These "Neganov-Luke" phonons [22] contribute to the 
total observed phonon signal, Ep, yielding 

Ep = Eii + eVi,NQ=ER+^EQ, (7) 

where Vb = 3.0 V is the bias voltage across the detector. Since Eq = En for electron recoils with full charge collection, 
Ep = (l + sYb) Er for these events. 

In this, as in previous low-energy analyses |121 I27j. the recoil energy for each event is determined from the phonon 
signal alone, after accounting for the Neganov-Luke phonons using Eq. [Tj This recoil energy estimate avoids incorpo- 
rating the poorer signal-to-noise of the charge measurement at low energy, but requires an assumption of the recoil 
type to correctly estimate the recoil energy. 

For electron recoils, the phonon-based recoil energy scale, Eer{Ep), is directly calibrated in situ using the 1.3 keVee 
and 10.4keVco Ge activation lines visible after neutron calibrations, ensuring that this energy scale cannot be overes- 
timated at the 90% C.L. Due to the calibration procedure, this scale is given by definition as 

Eeb.{Ep) = ^, (8) 

where we have used that eVb/e — 1 when biased at the calibrated voltage H = 3.0 V. 

The energy scale for nuclear recoils, Ej^r^Er), reported in units of "keVnr," uses the same calibration of the phonon 
energy scale as for electron recoils, but with a smaller correction for the Neganov-Luke phonons produced by a nuclear 
recoil. The mean ionization energy for nuclear recoils Eq,nr{Ep) is determined over the relevant energy range from 
the distribution of nuclear recoils from in situ ^^^Cf calibration data. The nuclear recoil energy scale is then 

Enr{Ep) = Ep - Eq^nr{Ep). (9) 

As stated above, the electron-recoil energy scale Eer{Ep) is directly calibrated with lines of known energy. The 
nuclear recoil energy scale cannot be directly calibrated in a similar fashion due to the lack of spectral features in the 
energy range of interest. Errors in the nuclear- recoil energy scale could be introduced by errors in the Neganov-Luke 
phonon correction relative to electron recoils, or dependence of the athermal phonon collection on recoil type. Since 
the Neganov-Luke phonon contribution accounts for only 20% of the total phonon signal and the ionization produced 
by nuclear recoils is well-measured using the ^^^Cf calibration data, the correction for the Neganov-Luke phonons 



does not contribute significant uncertainty to the recoil energy scale [T^]. The absolute nuclear recoil energy scale is 
constrained by comparison of the measured ionization yield in CDMS to previous measurements of the ionization yield 
for nuclear recoils of known recoil energy (see e.g. [28] and references therein). The measured yields are inconsistent 
with an overestimate of the nuclear recoil energy scale in CDMS, assuming the same ionization collection efficiency 
as previous measurements |12j . 

CoGeNT must similarly assume either a nuclear-recoil energy scale or an electron-recoil energy scale. In this 
analysis, we follow the prescription |3] of the CoGeNT collaboration to relate the ionization energy (calibrated using 
electron-recoils) to the nuclear recoil energy 



^Q = (0.19935)^] 



(1.1204) 
NR ' 



(10) 



good over the energy range 0.2 keVnr < i?NR ^10 keVn 



The CDMS II Exposure 

As noted in the main text, the CDMS WIMP-search data are not continuous over the nearly two years of exposure 
considered here, but include gaps due to neutron calibrations, warming of the detectors, and data periods removed 
due to data-quality diagnostics. Since data-quality diagnostics remove individual detectors, different detectors serve 
as WIMP-search detectors for different times. Figure [5] shows the history of exposure for the period considered on 
a detector-by-detector basis. The detectors were arranged in five "towers," and are identified by their tower number 
(T1-T5) and by their (top-to-bottom) ordering within the tower (Z1-Z6). 
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FIG. 5. The CDMS II exposure, displaying the detectors and time bins used in this analysis. For each time bin, a detector is 
colored blue if this detector's data was used in this analysis. Divisions between "runs" represent at least partial warm-ups of 
the dilution refrigerator used to cool the detectors. In order to avoid the effects of Ge activation, 20-day periods were omitted 
(red) following each ^^^Cf calibration time. 



Event Selection 
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FIG. 6. Event selection for the eight detectors employed by this analysis, shown in an ionization-energy vs. phonon-energy 
plane. The energy range of this analysis is indicated by thick green lines, along with borders between smaller energy bins. 
In this and following plots, the 5-11.9keVnr range has been subdivided into three equal parts: [5-7.3keVnr], [7.3-9.6 keVnr], 
and [9.6-11.9 keVnr]. AH events shown have passed data-quality cuts, as well as the veto cut and the Q-inner cut. Crosses 
represent events registering in multiple detectors ("multiples"), filled markers represent events registering in only one detector 
("singles"), and colored markers represent singles lying within the ±2a" nuclear-recoil bands, defined through ^^^Cf calibration 
independently for each detector and run. The edges of these run-by-run nuclear-recoil band definitions are also indicated. 



Limits on Modulation of Efficiencies 
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FIG. 7. Confidence limits on the relative amplitude and phase of annual modulation in the nuclear-recoil band cut efficiency 
(black) and the Q-inner cut efficiency (grey). Three different energy bins are shown, along with the total energy range (lower 
right). Contours are 68, 95, and 99%. 
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Rates vs. time 
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FIG. 8. Residual event rate as a function of time, for three event populations: multiples (top), singles (middle), and singles in 
the nuclear-recoil band (bottom), as defined in the text. Two energy ranges are shown, in the left and right columns. Because 
the multiples and singles populations are dominated by electron recoils, an electron-recoil energy scale has been used for these 
rates. 
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Modulation Spectra: Nuclear-Recoil Singles 
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FIG. 9. Amplitude of modulation vs. energy, showing maximum-likelihood fits where the phase has been fixed and the 
modulated rates M have been determined for both CoGeNT (light orange circles, vertical bars denoting the 68% confidence 
intervals) and CDMS (dark blue rectangles, with vertical height denoting the 68% confidence intervals). The phase that best 
fits CoGeNT (f06 days) over the full CoGeNT energy range is shown on the left; the phase expected from interactions with a 
generic WIMP halo (152.5 days) is shown on the right. The upper horizontal scales show the electron-recoil-equivalent energy 
scale for CoGeNT events. The 5-11.9 keVnr energy range over which this analysis overlaps with the low-energy channel of 
CoGeNT has been divided into 3 equal-sized bins (CDMS) and 6 equal-sized bins (CoGeNT). In the right plot, we also show 
the DAM A modulation spectrum (small grey circles), following the method of Fox et al. [29], for which we must assume both 
a WIMP mass (here, mx=10 GeV/c^) and a Na quenching factor (here, gjva = 0.3). Lower WIMP masses or higher quenching 
factors can push the DAMA modulated spectrum towards significantly lower energies. No attempt has been made to adjust 
for varying energy resolutions between the experiments. 



Rate Modulation: Nuclear Recoil Singles 
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FIG. 10. Feldman-Cousins allowed regions in a polar projection of modulated rate M vs. (j) for CDMS singles passing the 
nuclear-recoil cut (dark blue) and for the CoGeNT data (light orange). Three different energy bins are shown, along with the 
total energy range (lower right). Contours are 68, 95, and 99%. 



Rate Modulation: Multiples and Singles 
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FIG. 11. Feldman-Cousins allowed regions in a polar projection of modulated rate M vs. (j) for two event populations dominated 
by electron recoils: multiples (light blue) and singles (dark red). Three different energy bins are shown, along with the total 
energy range (lower right). Contours are 68, 95, and 99%. 



